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PRACTICAL CALCULATION OF LAMINAR AND TURBULENT 
BLED-OFF BOUNDARY LAYERSl 

R. Eppler 

Bolkow-Entwicklungen KG 


1. Introduction 


/221* 


Bleed-off of boundary layer naterial through snail ori- 
fices or porosities at the wall in contact with flow has long 
been known as effective neans for inproveaent of flow character- 
istics [1]. On the one hand it pemits longer conservation of 
the laninar boundary layer and on the other it helps prevent 
separation of the turbulent boundary layer* or retards it (2], 
both of which can lead to considerable reductions of drag. 

Often the two effects of bleed-off overlap. Where bleed off 
is present in an air foil, for preservation of laninar ity at 
high-speed flight, it will also have a certain effect in re- 
tarding the separation when it can no longer prevent the 
reversal during slow flight. Conversely, continuous siphoning 
for lift increase, i.e., for retardation of separation during 
slow flight, must in principle have a favorable influence on 
the reversal at high-speed flight during a change in pressure 
distribution. Up til now hardly any thought has been gi^'en to 
the combination of these two effects, since entirely different 
amounts of bleed-off are required at different points. But 
lately there has been some approach to this combination. A 
valuable aid in realizing such a bleed-off combination is a 
method of calculation that permits uniform treatment of all 
cases . 


^Report from Bolkow-Dcvelopments, Inc. 

*Numbers in the margin indicate pagination in the foreign text. 
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This study will report on a solution of the probleu. In 
the foreground will be the uethod of calculation with one paraa" 
eter by aeans of Kimdn*s [3] Conservation of Mouentuu and of 
Ifieghardt's [4] energy theoreu, which was first applied suc- 
cessfully by A. Valz for lauinar boundary layers [S]. Since 
then the uethod has proven itself for uany other purposes. 
Wieghardt [6] calculated lauinar boundary layers with bleed- off 
in which the uininuu auount of suction needed for conservation 
of lauinar flow was also calculated. Truckenbrodt [7] extended 
Walz*s uethod to turbulent boundary layers by using seuieupir- 
ical equations based on experinental and theoretical investi- 
gations by Ludwieg and Tilluann (8]. and Rotta [9], instead of 
different relations derived frou velocity profiles during 
lauinar flow. The uatheuatical siuplifications used did 
result in unattractive errors for a few special cases , partic- 
ularly in the vicinity of the turbulent separation point. 

Schulz [10] offered an iuproveuent, particularly effective for 
that point. A procedure by Walz []]] for the turbulent case 
including coupressibility, which operates without mathematical 
siuplifications, is seuigraphical and thus little suited for 
the digital computer. Only lately have Schlichting and Pechau 
[12,13] reported on the extension of Truckenbrodt *s theory of 
turbulence to bleed -off, in which additional mathematical 
simplifications were carried out. Such simplifications are 
particularly annoying since the empirical fundamentals of the 
theory of turbulence are still uncertain. It is hard to dis- 
tinguish whether existing shortcomings are of purely mathe- 
matical nature or if they are due to the empirical fundamentals. 

For a case without bleed-off it was recently tried to 
avoid the simplifications introduced by Truckenbrodt; this led, 
simultaneously, to further unification of laminar and turbulent 
cases and a possibility for simple programming on digital com- 
puters [14]. It was demonstrated that the empirical 
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fundaaentals used by Truckenbrodt are better than the results 
gained with the Method had led one to suspect. Extension of 
those nethods to the case of bleed- off, which will be presented 
subsequently, again pemits sinultaneous treataent of laninar 
and turbulent boundary layers and solves the approxination dif- 
ferential equations, according to the best presently available 
enpirical fundaaentals, without an uncertain error. Here, too, 
a big difference is shown fron Schlichting and Pechau. Unfor- 
tunately, possibilities for conparison of experiments are still 
sparse. tut the results gained are entirely plausible and 
because a large amount can be obtained in a simple way, and 
quickly, it has become possible to check any and all experi- 
mental results with them. Differences between theory and ex- 
periment can be pursued back to the empirical fundamentals 
because of the easily grasped method of solution, permitting 
their verification and correction. 

The case without bleed-off, which is completely contained 
as special case, will be treated in detail once more, since my 
cited lecture appeared in print without the corresponding illus 
t rat ions. 

2. Basic Equations 


The following quite well known relations (with subse- 
quently explained symbols) can be formed from L. Prandtl's 
well known partial differential equations for velocity u(x,y) 
in the boundary layer (see Fig. 1), by formation of moments 
through application of Bernoulli's equation for the outside 
flow U(x) in the plane incompressible case, according to von 
Karmdn and Wieghardt: 
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1 m 


(1) 



Fig. 1 Sketch of Coordinates 
and Velocities. 

piacesent thickness is shown as 


!* *00 + u jx 

I a MwnjL\ - 4 _ r* 


( 2 ) 


where is non»l velocity 
at the wall, positive in the 
y-direction (blow-off), . 
According to Walt [S] dis- 



(3) 


the so-called pulse and energy loss thickness as 


and 



(4) 


( 5 ) 


Symbol to, which occurs in the Conservation of Momentum equa- 
tion (1), is shearing stress at the wall 


in 


( 6 ) 



(7) 
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fro« equation (2), generated at first formally during formation /223 
of moments » in which d is interpreted as energy dissipation and 
t as change in turbulence (disappearing for laminar flow). The 
whole equation (2) can be treated as energy theorem. The (con- 
stant) density of the flowing medium is f , kinematic viscosity 
is V. 


Shearing stress t from equation (7) is analogous to (6) 
for laminar flow 


bat still represents a big uncertainty for the turbulent boun- 
dary layer. 

Equations (1) and (2) are at first valid only for the as- 
sumption that the solution u(x,y) of the partial Prandtl dif- 
ferential equation is known for the given functions of U(x) 
and v»(x) and that the values required for equations (1) and 
(2) can be calculated from (3) to (7) with it. For the case 
of turbulent flow velocities averaged over time can be used 
during consideration of the contribution of turbulence to 
shearing stress. 

Reversal of the calculation procedure for an approximate 
solution is carried out by making certain assumptions about the 
distribution u(x,y), and for the case of turbulent flow also 
..bout shearing stress x. This will be described in more detail 
in the following section. 

3. Form Parameter, Wall Shearing Stress and Dissipation for 

Turbulent Plow 


Two relations, (1) and (2), are available for A 

statement can be made for u(x,y), which contains two independent 
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values. Usually one independent value is taken for the thick- 
ness of the boundary layer, that is for distortion of u(x,y) in 
the y-direction, another for the shape of the boundary layer 
through a stateaent like 



( 9 ) 


by entering this arrangeaent into (3) to ( 8 ) one gets all the 
values defined by (3) to (7) as functions of fi(x) and X(x); 
equations ( 1 ) and ( 2 ) becoae a system of differential equations 
for these two unknown values. Approx iaat ions can be obtained 
subsequently with their solutions, for u(x,y) from equation 
(9), for the 6 . from (3) to (5). 

There is a large number of different equations (9) and 
many methods and procedures for their solutions. The following 
one, introduced by Walz and Wieghardt, appears to be the most 
easily understood calculation method. When entering (9) in (3) 
the transformation n*y /6 shows that 6 i is proportional to 6 , 
i.e., that 6i/6 depends on X only because 


^1 = 

9 


( 10 ) 


The same holds for and 63 . This makes the quotient of two 
5^ a function of X only, which can be obtained once and for all 
for each given (9) . For instance, by putting 

I7 ■ (11) 


and 


61 
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( 12 ) 



Hi 2 *F(Hsz) can be obtained by elimination of A and 6i, which is 
needed in equation (1), expressed in 6t and 6$. Should one 
succeed in expressing shearing stress and dissipation similarly 
simply through 62 and 6s, the solution of the (1) and (2) 
equation system is recommended without further transformation. 
According to (6) and (9) shearing stress at the wall is 


u_ 

p If* 


r(f\ 


I A. 




(13) 


Here is 62 / 6 -, already shown as function of A only; this can 
also be arranged for 3f/3 n through (9). Through reversal of 
( 11 ) shearing stress at the wall can also be expressed through 
6 2 and 63 . 

Dissipation can also be brought into suitable form through /224 


<1 t I 


oo 



0 


(14) 


As soon as the functions Hs 2 (A), HisCHss), ^^(Has) and D*(Ha 2 ) 
have been obtained for a given equation (9) the system (1) and 
(2) must be integrated numerically. With the results 62 (x) 
and 5g(x) one can find Haz, from it A and, with the aid of 
equation (9), the basic velocity distribution u(x,y). In most 
cases the last step is not carried out; generally it is quite 
satisfactory to have in 62 a value for friction losses and in 
A [or equivalent value, like Haz froir. (11)] an indication for 
the shape, resp. the separation tendency of the boundary layer. 

Everything has not been traced back to equation (f*) . The 
closer on. comes to the exact solution u(x,y) with it the better 
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will be the results. The aost iaportant boundary conditions 
are, therefore, fulfilled in advance in function f(X,y/6) froa 
equation (9) 

f(X,0) - 0, f(A,«) - 1, -4 (15) 

y 

It is already sufficient to fulfill the last two boundary 
conditions in the finite, i.e., for the distance of a certain 
boundary layer thickness 6 froa the wall, or n»l. Even though 
Walz has shown that after consideration of (15) the special 
equation (9) has little influence on results where no bleed- off 
occurs, soae attention will still be devoted to that question 
since the joint treataent of bled off boundary layers and unin- 
fluenced ones includes a large range of various velocity pro- 
files, which could create inaccuracies. Important is only the 
influence on the three functions required for (1) and (2). 

A particularly good premise was already provided by 
Pohlhausen (15) . He puts 

S • Pi(n) + A(x) P 2 (n) for n » 1, 

g = 1 for n > 1 

and fulfills the boundary conditions with the fourth degree 
polynomials Pi and P 2 (equation (15)). This so-called Pi,- 
statement is, like the improved Pg -statement , too coarse for 
our purposes. Much better suitable are surely the velocity 
profiles that are generated at the so-called "similar boundary 
layers" (Hartree profile), particularly for flows where the 
form parameter does not vary quickly. Walz [5] already em- 
ployed these solutions, which are exact for certain velocity 
distributions, for calculation of the required connections. 


8 





Inaccuracies that have slipped in do not play an important 
role in Walz's case but they carry more weight in the calcula- 
tion of bleed-off conditions. For that reason a few more 
similar solutions were obtained with the corresponding values 
for Hs 2 » Hi 2 » and (Hartree profile). In Fig. 2 data 
points are plotted for the similar solutions without bleed-off 
for increasing pressure (Hartree profile) and others for solu- 
tions of the flat plate with bleed off. In addition a data 
point is plotted for H 3 2 * 5 / 3 , which was obtained from the so- 
called asymptotic solution of the plate with constant bleed- 
off. The following approximation was chosen for the calcu- 
lation to provide good approximation for the entire range: 


a. Values obtained from the Hartree profiles were approx- 
imated for the range 1.51509»H32*1. 5/258, which corresponds to 
the pressure increase range between the flat plate and separa- 
tion for Hartree profiles. Hsa*!. 51509 is, therefore, the 
separation boundary. 

b. In place of the Hartree profiles of pressure decrease 
the similar bleed-off profiles of the plate were used, comple- 
mented by the point of asymptotic bleed off. 

Boundary layers free of bleed-off are characterized somewhat 
less accurately that way. But in this case accuracy is of less 
importance. 

The approximations 


//„ = 4,02922 - (r,8.1.60l«2 - 72I..'.39I6 | 227,l«220 ^ //^ |.r,ir.09 

f»r i,r.ir>09 //„ i,r,72.',« , 

//„ =r 79,87081.', — 89, .'■,82 142 //„ ^ 2.'i,7l 5786 //*, 

fOr //„ > l,.57258 , 


(16) 


/225 
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( 17 ) 


e* =s 2.S12589 — 1.68609S f/„ + 0,391S41 llfj — 0,031720 /|J, 

for 1,51509 1,57258, 

e* = 1,372391 —4,226253 //« + 2,221687 If^ 

for /fa, >1,57258, 


If* = 7,853976 — 10,260551 //„ + 3,418898 1/^ 


( 18 ) 





Fig. 2 The General Boundary Layer -Functions Hi2(H32), e*(H 32) 
and D*(H32) in Several Directions. 
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OP POOR QUAIOT 


are suitable for generation of computer programs. The dif- 
ference between these approximation functions and the approxi- 
mated values is smaller by an order of magnitude than the dif- 
ference between the various boundary layer profiles being dis- 
cussed. For calculations by hand it best to use accurate 
diagrams of the functions (16) to (I j, which are easily ob- 
tained with Table I. 

The choice made provides very good coverage of similar 
solutions without bleed- off with pressure increase to the flat 
plate (Blasius Flow ) , as well as for similar bleed-off boundary 
layers of the flat plate. All other calculations of boundary 
layers with the present method contain an error brought about 
purely by the procedure, which will be explained with several 
examples . 


Table I 


Values of Laminar Boundary Layer Functions 
Hi 2 (Hs 2 ), e*(H32) and D*(H32) 


«.l 


i 

li* 

I.S15(»9 

4.02922 

j 0,018100 

0,156.38 

1,51576 : 

3,83366 

; 0,01107 ! 

0,15639 

1,52099 

3.48079 

’ 0.049 U 

0,13672 

1,52630 

3,29673 

0,07303 

0,15711 


3,02095 

0,11796 

0,16006 

1,54803 

2,87178 

0.14828 

0.16292 

] .55568 

2,77120 

0,17166 

0,16371 

1.56214 

2.69668 

0.19071 

O.ltiBll 

1, 56771 

2,638.38 

0,20674 

0.17100 

1,57258 

2.59110 

0,22052 

0.17310 

1.60353 

2.3(611 

(1.30706 

0.19219 

1,62256 

2.224^9 

0.36108 

0.20.08 

1,6(289 

III 182 

(U26H5 

0.22539 

1.66020 

2.02B46 

0,18039 

0.21313 

J, 66667 

2.00000 

il.SOOOO 

0.25000 
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Fon Parmieters. Mall SLearing Stress and Dissipation for 
TfirKUTentTloif : 


Since the laninar boundary layer is easily handled, ac- 
cordinf to the results in section 3, % corresponding path will 
be followed for solutions of the turbulent layer. The only 
significant difference will be that shearing stress can no 
longer be expressed through equation ( 8 ) in the turbulent 
boundary layer. As function of the turbulent notions of vari* 
ation and the energy transport connected with it additional 
**apparent** shearing stresses occur at first, which call for 
other relations in place of ( 8 ). In addition, the velocity 
profiles take on conpletely different foms so that statenent 
(9) does no longer cover then conpletely and energy dissipation 
cannot be obtained the sane way as before. But wall shearing 
stress ( 6 ) is not influenced since no turbulence is possible 
and effective in close proxinicy to the wall. 


So far it has not been possible to get nseful results for 
the theory of turbulent boundary layers based only on laws 
about the turbulent shearirv stress together with velocity 
statenents. But just as we could liberate ourselves, with the 
aid of general iunctiens fron the velocity profile, of the boun- 
dary layer in the laninar case so will relations which pernit 
calculation of wall shearing stress and energy dissipation fron 
the independent variables ,62 and 6 ], also suffice in the turbu- 
lent case. Such relations for the wall shearing stress were 
stated by Ludwieg and Tillnann ( 8 ] as 


»• 

•• 


0.123 
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(19) 


for d issipation by Rotta (9]. Rotta's results were used by 
Truckenbrodt [7] in simplified form 
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4 + 1 


( 20 ) 


‘ MtUj 




JKicli h«s proven itself well for the case without bleed* off ( 1 ) 
like the connection of fom paraneters based on Wieghardt (16] 


( 21 ) 

The validity range for equations (19) to (21) nust be ex- /227 
plained in nore detail since for the case of bleed- off larger 
fom parameter ranges will again overlap. At first (19) was 
obtained for values of Hia between 1.2 and 2.4, corresponding 
to H |2 values between l.S and 1.85 per equation (21). Equa- 
tion (20) also is valid in the sane range. No overlap in 
range will be of concern for the pressure increase up to sep- 
aration, which occurs at about Hs 2 * 1 . 5 . But a strong tendency 
to high Hs 2 values is present for the range of pressure de- 
crease or for bleed- off. Pechau [3] has given examples where 
Hs 2 values from 3 to S were reached. The range of high H 32 
values or low Hi 2 values nust, therefore, be observed sepa- 
rately. 

A possibility for extrapolation fron the Ludwieg-Tillmann 
law is offered by the following boundary consideration. After 
62 and 63 are defined by equations (4) and (5) one gets 








( 22 ) 
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as long as u<U» since the second parenthesis can be naxinized 
with the value of 2.« As long as there are no excess velocities 
present in the boundary layer, which is not to be expected 
without tangential blo%r-off, 

H,a < 2 (23) 

according to statenent (11) . Based on the sane assumption 
Hii>l fron comparison of (3) and (4). The nature of boundary 
layer profiles leading to Hsi values near 2 is also easily 
grasped. For instance, if one wants to reach 

Hj2 ^ 2 * € (24) 

then for l-u/U*u in statement (5), 

00 

6 , • |(l-u)u(2-u)dy > ( 2 -e )62 (25) 

SO 

oo ® 

ju*(l-u)dy < c 62 • e|u(l-u)dy. (26) 

As long as u> 0 , i.e., as long 
as no excess velocity is 
present, this would be valid 
if u<r- everywhere. This 
surely is not true close to 
the wall though. If we as- 
sume, according to Fig. 3, 

that in a range close to the 

Fig. 3 Sketch for the Esti- where 0*y»yx a share of 

mates Carried Out. ' ' 

621 is already present of 

pulse loss thickness while u*e in the rest of the range then 



14 



(27) 


5* ■ ^2 1 ♦ F 

when F, the share of displacenent thickness 6 i,is outside of yi. 
The left side of (26) can be estinated at the lower Unit as 

/■»(! — 5)iir > j«*(i — ••) *y “ /‘(s — «)5(i — +e/'s (1 — S)^ = (JC + e)^ , (28) 

Q ft O ft 

where k is probably a constant different fron 0 and snaller than 
1. Using estimates (27) and (28), (26) becones 


(29) 


or 


{K + t)d„<e(6n+ Fi 

F>1 


(30) 


Whenever a finite 621 and a finite ui is present for any , 
Hs 2 can only go toward the value of 2 when i.e., when the 

share of displacenent thickness far away from the v n, where 
u<e, increases toward infinity. Simultaneously witi. ( 2 -H$ 2 ) , 
(Hi 2 - 1 ) also decreases towards zero because 



’ Jl 

J • 

0 

J 

0 




(31) 


Here I /62 decreases to 0 at least with e for 2-H32*e-*-0, be- 
cause of (30), while F /62 is limited. The entire right side 
of (31) decreases toward 0 with c. From this follows 

lim Hi 2 ■ 1 (32) 

H32 -»• 2 


original 
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This relation is not considered in the interrelations used so 
far (compare (21) )» but can be introduced after some snail 
ch.nges that are in the uncertainty range of the entire inter - 
r<> lation. For bleed-off calculations 


Hi2 


11 H|2+15 
48 Hs2'35 


is recoanended in place of statenent (21). 


(33) 



Fig. if VariiUj Interrelations 
of Form Parai.. ters Hi 2 (H 32 ) 
for a Turh'^lent Boundary Layer. 


the . ncertainty region. 


According to Fig. 4 there are 
hardly any differences between 
(21) and (33). They are ef- 
fective only in the vicinity 
of Hs 2 * 2 » in connection with 
the following derived basic 
laws for wall shearing stress 
Other form parameter func- 
tions are also plotted in 
Fig. 4 for the sake of clar- 
ity. Significant differences 
are present only where H 32 
values are low. While the 
interrelation used by 
Truckenbrodt [7] provides 
slightly larger values of 
Hi 2 , lower values are gener- 
ated by another interrelation 
based on data by Cornish, for 


\fter extensive calculations with all three essentially 
different interrelations, which were partially reproducevi by 
Seyb fl7j, the interrelation (33) was recognized as being use- 
ful In the normal range as well. 
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Theorea (19) for the wall shearing stress can also not be 
extrapolated in the present fora for bleed-off exaaples. Ac- 
cording to (L') wall shearing stress decreases toward zero for 
Htx*^l because of (30), since 6x increases toward infinite. 

This aust be avoided. Apparently, very large contributions of 
F to fit or 62 aust be available for saall values of Hix-*-!, de- 
rived frcMi the regions far away froa the wall with very saall 
differences in potential flow. Those regions cannot have any 
influence on the wall shearing stress. A tern aust be added to 
(19), therefore, to take care that the values for 2-H32*e, 
which increase as 1/e, not have any influence on 62. That is 
done through 


The nuaerical constants in this relation are fixed so that 
there is the best possible agreement with the wall shearing 
stress theorems used so far, in the region of the most fre* 
quently occurring values of Hi 2. In Fig. 5 the comparison with 
(19) is shown for three values of U62/V vs. Hi 2. The difference 
is appreciable only for values of Hi 2 close to 1. For that 
region (19) was not derived. The boundary condition treated 
in (34) is represented here. 

A further comparison of the wall shearing stress theorems 
is shown in Fig. 6. In addition to (19) lines according to 
Rotta [9] and Cornish [18] are plotted. The new theorem (34) 
is not essentially different from the others in the normal 
region. The straight lines are a little less inclined compared 
to (19)3 which is caused by the slightly lower exponent of the 
Re number. That appeared expedient so as to improve the agree- 
ment with Rotta and Cornish for higher Re numbers. 


/229 
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Equation (2) for dissipation has been very reliable. It 
can already be seen for the laminar case that the fom paran- 
eter has less influence on than on e*. In the turbulent 
case it has so far been completely eliminated (according to 
( 20 ). 



Fig. 5 Turbulent Wall Shearing Stress as a Function of Hi 2 . 



Fig. 6 Turbulent Wall Shearing Stress as a Function of Re^^ 
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This can only be naintained approxiaately in the following 
treatnent. since dissipation for H3 2'*^2 or Hi2->1 aust also only 
be dependent on (U62/v)^(Hi2*l)« which is checked by aeans of 
(7) . Instead of (20) 

“U ‘ ^ ‘ (35) 

will be used subsequently, where the boundary condition is 
fulfilled and only an insignificant influence of Hi 2 is present 
in the nomal region. A large nuaber of coaparative calcula- 
tion bas furthermore established that dissipation has less 
influence on the results than shearing stress and that the 
boundary condition is aore important than a really exact value. 
For all normal values of Hi 2 (35) is, therefore, equivalent to 
( 20 ). 


5. The Numerical Treatment 


/230 


The differential equations (1) and (2) are completely de- 
fined for the laminar and turbulent case, with the relations 
derived in sections 3 and 4, and need only be integrated numer- 
ically. This can be done with known methods. Wherever the 
integration process will be discussed subsequently, it will be 
only to point out various characteristics of the differential 
equations themselves, which are used advantageously for the 
numerical treatment. Solved according to the derivations of 
the desired functions, equations (1) and (2) are as follows 


dx 

dx 


= ~(2 I //„) 


#>2 dC 
V dx 



3 A 3 fIV I t^d ^ t r, 
II dx ^ ~ ' V 


-}■ 




(36) 


( 37 ) 
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eV lUa, 

(buntiiar)* 

(38) 

= 0.04S716((Jftt— 1) «» 

(lurbttleni) • 

(») 

*«!/»- JK 

(lainiaar) » 

(40) 

2 7 ^ = 0.0100((l/„ - 1) • 

(turbuiriit) . 

(«) 


As abbreviation the Reynolds number formed with 02 and U 


Re 



(42) 


was introduced. The function Hi 2 needed in (37) and (38), 
is defined by (16), resp. (21), e* and D* by (17) and (18), 
resp. Table I. 


Equations (36) and (37) are apparently dimensionless since 
velocities or lengths appear only simultaneously in numerator 
and denominator. It is all the same, therefore, in which units 
velocities and lengths are being measured. Vie introduce as 
unit for the lengths 1*1 and for the velocities U^® 1 » both 
chosen to be most suitable for the particular task. All lengths 
and velocities can consequently be considered as ratios of the 
chosen units, resp. as functions made dimensionless along with 
the units. This changes nothing in the entire formula system. 
Only in the introduction, for convenience's sake, of the Re 
number formed with the units 


Re ■ 


1*1 M 

V V 


(43) 


will in place of (42) 


Re 


62 


Re U 62 (»— 


U6, 


( 44 ) 


be used in the future. 
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The tt^fferontial equations (36) and (37) have the forn of 
• ^(x»52»6s)» ■ 8 (x» 62 ds)» 

where the only difference between laninar functions f and g in 
the laninar and turbulent cases is that for wall shearing 
stresses, energy dissipation and form parameter^ sometimes equa- 
tions (38), (40) and (15), and sometimes (39), (41) and (21), 
must be used. U(x) and Vo(x) are given functions, dU/dx can 
easily be formed during integration. 

For step-by-step integration Runge's second order method 
is most suitable. With initial values Xi, 621, 892, Ui, vqi, 
final values X 2 >xi<*-Ax, and analog U 2 *,v «2 given, the differences 
between the desired variables are calculated by means of 

Ax , = <?(*!• 4|»’4ji) -1* 

f 5 I*. *- = 

For hand calculation little schedules for the calculation /231 
of the right-hand sides f and g must be made. Often it will be 
noticed that after the first half-step a repetition of the cal- 
culation will have no significant influence and so this step 
can be eliminated. 

On occasion one will observe another, much more unpleasant^ 
characteristic of the solution. Differences fluctuate in a 
strange way rather heavily around the true value, particularly 
when the boundary layers are still thin, as for stagnation 
point flows or in the vicinity of an intake. The reason for 
it is that neighboring solutions, to which one turns because of 
inaccuracies or mistakes in rounding off, tend strongly toward 
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the desired solution, particularly for H 92 * The step>by-step 
integration follows, according to Fig. 7, by going halfway in 
the direction toward the tangent to the point of departure 
then calculating the tangent direction and following in the 
tangent direction for completion of the entire step. When the 
neighboring solutions tend too strongly toward the solution, as 
sketched in Fig. 7, the error during one step increases greatly. 



Fig. 7 Sketch of the Numerical 
Instability Occurring. 

For that reason alone it was not 
studies . 


The reason for niu&erical 
instability is not based on 
this particular approach to 
a solution but is character- 
istic for differential equa- 
tions (36) and (37) and, 
therefore, in more or less 
apparent fashion built into 
all previous studies based on 
the pulse and energy theorem. 
The better the applied numer- 
ical method, the more the 
instability becomes apparent, 
discovered in the previous 


It is apparently of little use to employ a method of a 
higher order, like that of the fourth degree by Runge-Kutta, to 
remove the instability. In that case the instability can even 
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greatly increase. The successive step decrease has so far been 
the Bost reliable, generally applicab leeway. As long as the 
nuaerical calculation is carried out by hand and one can survey 
not only the fluctuations but the entire calculation' constantly 
other paths are open. From time to time a check can be made 
whether the initial calculation, for stagnation point or intake, 
for instance, as discussed later, can be further extended. 

Often that can already lead away from the instability region 
but it must not be considered as panacea. An uncritical ex- 
tension of the initial step can lead to serious errors. When 
programming instructions for calculations on digital computers, 
which are very suitable for it, checks must be built into 
every step, for recognition of numerical instability. When 
it occurs only a decrease of the step length is still indi- 
cated, best to half. The required criteria will now be pro- 
vided, briefly. All other details of programming need no 
further explanation. We designate a few of the functions that 
occur in the course of the integration step 


Ha 2 i 




Haaz 


^ S 2 

522 


After that there should be stabilization with 


«) 

W 

r) 

d) 


— 2 lljg -f J/„,| > f| (f, » 0,0001 bis 0,01) , 
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The excluded possibilities fron (b) to (e) occur nornally any- 
way only in cases that already fulfill the condition (a) and 
thus aim for a step decrease. They must be recognized before 
check (a) only because they might require [as in (39)] useless 
calculations which would lead to a stop of the computer before 
check (a) . Once the step decrease has been programmed 5 " s 
no trouble to catch all these absurd cases and to chan^ 'r 
to step decrease if need be. 

Quite a few other refinements '‘f the calculation can be 
achieved through step decrease. For instance, Ha a often de- 
creases steeply in the vicinity of separation points, creating 
some inaccuracies with it. This makes it advisable to limit 
the entire contribution to change of H 32 ^nd to carry out the 
step bisection when 

f) is |H322-H32 i| > 62 (£2 * 0.01 O 0.03). 

Naturally, the bisectioned steps are checked again for 
fluctuations of differences that are too strong and bisectioned 
again^if needed. This bisectioning process shuuld be carried 
on at random, a lower limit being set for Ax»(1/2)*'Ax which, 
in order of magnitude, is I/IOOO of the smallest initially oc- 
curring step size. Going below this limit may stop the program. 
It has been shown that a large number of programming errors can 
be laid to continued step bisectioning. These errors are found 
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easiest after the program is stopped because the step 
have become too small. If one were to try to continue r/ne cal> 
culation the search for errors would at best be made more dif- 
ficult, since memory locations are overwritten whe: e errors 
could otherwise be recognizer. 

*'*uring bisectioning of the examples shown later U and vo 
were lijiearily interpolated, i.e., held constant over the en- 
tire initial step dU/dx, even though it was subdivided into 
perhaps 32 or 64 partial steps ~Ex. This is justifiable when 
the step size of the particular problem was chosen so small 
that linear interpolation between supporting points is mean- 
ingful within the context of this problem. With the integra- 
tion method chosen it leads to a second order approximation 
when the output interval Ax is handled without bisectioning. 

A step decrease is, after all, introduced for elimination of 
purely numerical instability and not for increase of calcula- 
tion accuracy. It is easiest to increase the number of support 
points if there is concern about the linear interpolation. 
Better interpolation of U and dU/dx is not a principal diffi- 
culty, however. It is to be used with care only when dU/dx 
shows nonuniformities, which happens often in practice. 

The numerical integration of differential equations has 
now been taken care of. Calculation by hand is very simple; 
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just the saue it is aore effective to prograa it so ms to ob- 
tmin fmr more voliminoos result aateriml mt considerably lower 
cost. A few natters have to be considered, as always during 
prograaning, which do not becoae conspicuous during calculation 
by hand since there is constant monitoring of the calculation 
flow and innediate intercession when something is not quite 
normal. Now there remain only a few special cases^at the start 
of the calculation, at the boundary layer shift and at the 
separation, t'> be discussed. 

6. Initial Values 


The integration of the differential equation system (36), 
(37) , so far discussed, does not assume disappearing values for 
U, <x and 6 %. These conditions are not fulfilled at the start 
of a boundary layer. The two most important cases are the 
**shockfree approach** of the potential flow at a^^vp edge and 
the stagnation point flow. In both cases cons iddfat ions of the 
boundary are required, the first steps of the boundary layer 

calculation being easily carried out with their results. 

$ 

a. Shockfree Approac h 

In this case the boundary flow starts with a U that does 
not disappear. It is shown through well known procedures that 
the boundary layer then starts with disappearing 5x and 6$. In 
(36) and (37) only those terms with the smallest powers of fix 
and 5) are of importance for small values of fix. If the others 
are neglected then 



•’ !»• . 


(45) 


(46) 
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BqiMtioB (46'i provides a value for Haa» which is exactly the 
oae that also holds true for the plate flow U-constant. Since 
we have obtained the interrelations of the fom parameters 
(16) and (18) fron the sinilar boundary layers, which contain 
the plate flow, (46) oust be fulfilled for the Has value of 
the plate flow, namely Hia*l.S7258. That can be checked easily 
with Table I and the equations (16) and (18). It also serves 
for easy calculation of 6t fron (45). By putting the value of 
e* that belongs to Hjs«l. 57258 into the equation at the onset 
of flow where x"0 and U«U«, 

• 0.66411’''^l2^ . (47) 

4$ is also known through (11) and the equations (36) and (37) 
can be further integrated numerically. 

b. Stagnation Point 

In the vicinity of the stagnation point U>U'x. It is 
shown that the functions 62 and 63 must be constant for small 
values of x and for v«-0. The well known boundary considera- 
tion then furnishes 


and 


(48) 

2(2 ; 

\ I 

(49) 


Although the boundary consideration proceeds independent of Vo 
during the inflow, voj^O would change equations (48) and (49). 
But since bleed-off is hardly used in the vicinity of the 


*From 48. See (17|. 
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stagnation point it is sufficient to neglect it for at least a 
small part of the way in initial calculations. 

If the same form parameter interrelations had been used 
for Hs2>1. 57258 as for the similar boundary layers free of 
bleed>off then (48) should be fulfilled for the exact value of 
Hj 2 at the stagnation point flow, namely Hs 2*1. '^2575, which 
would give 62«0.29234//RFTF. But since more consideration 
must be paid to bleed-off in this region, condition (48) is 
fulfilled for (16), (17) and (18) when Hs 2 «l. 61998, which leads 
to 


«2 - 0.29004 i (50) 

/Keir 

The deviation in 62 is unimportant, being less than It, 
but is more significant in Hsa since there is no way to compare 
the absolute amounts. When referred to the H 32 -region between 
separation (1.51509) and asymptotic bleed-off at the plate 
(1.66667) the deviation in H 32 amounts to 3.8t, which is toler- 
able in relation to other methods of approximation. 

This error is, incidentally, a basic one for all methods 
with one parameter. It stems from the fact that in reality 
different boundary layer profiles can belong to the same form 
parameter H 32 and so can different values for e* and D*, while 
in this case a special relation was selected. According to a 
suggestion by Seyb [1] one might attempt to use various func- 
tions Hi 2 » e* and D* for cases with or without bleed-off. 

That would certainly allow significantly impr.. i^ed accuracy in 
the handling of many other examples, for instance the stagna- 
tion point flow and also the flat plate with blow-off, which 
is particularly poorly covered in section 8 . It is question- 
able, however, if that will always involve higher accuracy, 
particularly for examples where there is only partial bleed-off. 
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The error at the stagnation point flow uy be taken as 
indication 6f the general shortconing of the Method, even though 
it is probleaatical to judge the efficiency of a Method froM 
special exaMples. ExaMples with Much SMaller errors will be 
calculated but soMe with considerably larger deviations can 
also be constructed. Such a one is the flat plate with blow- 
off. In practice it is uninportant but it clearly deMonst rates 
the liMits of the Method. 

7. Boundary Layer Separation, Boundary Layer Shift 

a. LaMinar Separation 

The separation point of the laminar boundary layer is de- 
fined by disappearing wall shearing stress. The disappearing 
wall shearing stress is assigned to the value H32«1.515P9, since 
a clear relation between the functions e* and Hs 2 was introduced 
in the present method. 

It is no problem for calculation by hand to fix the exact 
point where reaches the separation value, either by extrapo- 
lation or by change of step sizes. 

For computer programming there are again certain consider- 
ations, since H32<1. 51509 can never be entered into equations 
(16) to (18) because that would lead to - root with negative 
radicand and again to a computer stop. The following meihod 
has been shown to be the must reliable one. A check is made 
after each half-step whether H52’^l-S1509 or H 32 2<1. 51509. If 
yes, new values for Ax must be obtained by the Regula falsi 
(approximation method for obtaining the solution of an equation) 
with the aim to repeat the step with the changed i^x so that in 
the end H 322 *l • 51509+e. If the separation boundary is again 
transgressed aftei this step the Regula falsi is used again for 
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calculfttion of a.soaewhat saaller stO]> size, if Hst 2 « 1 . 51 S 09 
calculation is continued. It only renains to check before each 
partial step that H 9 z< 1 . 51 S 09 '*- 2 e. If so then the starting 
point of the partial step is the laminar separation point. 
Inunction e is only a threshold value, which prevents the pro- 
gram from getting lost in a rounding off error. At the same 
time the accuracy with whicn the separation point is determined 
can be gauged with e. With e*O.S . 10 '^ that accuxacy should be 
good enough. With this value very few step decreases are re- 
quired until the point of separation is found. 

If the calculation is not finished after the laminar sepa- 
rationtbut continued for the turbulent case, it must be 
remembered that the "Kx of the preceding calculation step did 
not result any longer from step bisect ioning. In contrast to 
step reduction through stabilization one would not come to the 
final point of the total step if one were to continue the calcu 
lation with the existing "Ex. For that reason the entire piece 
between the separation point and the end of the integration 
step is used as l^x after the laminar separation. After that, 
step bisectioning through stabilization is again permissible. 

b. Turbulent Separation 

Turbulent separation still presents a big uncertainty 
today. It is generally assumed that separation can occur for 
H| 2 * 1.58 and that it definitely has happened for H 32 = 1 . 46 . 

These boundaries can be considerably extended (upward and down- 
ward) particularly for bleed-off. Fortunately, the uncertain- 
ties are not quite as big as it may appear at first glance from 
the H32 values. According to Seyb [ 19 ], and to one of my inves 
tigations [20], the turbulent boundary layer does not permit a 
layer gradient any longer, if the value of H32 is permitted to 
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drop below 1.S8 (without bleed'off) . The result is» conversely, 
that for H)a<l.S8 the influence of the pressure gradient of the 
forn paraneter is very sensitive. In many cases the change of 
Ha 2 in that region is so great that only a snail distance lies 
between the two points at which the values of Hs2*1.58 end 
Ha2*l*46 are reached, there being little influence on the point 
of separation by the chosen value of the separation parnaeter. 
For bleed- off, however, the influence of the anount bled off 
is quite critical, so that one can stay on the safe side with- 
out mich additional bleed- off. 

All the exanples shown here were calculated up to the 
value of Ha 2*1. 46. A decision must be nade in each case 
whether another separation criterion would have had significant 
influence . 

The location of the point where Ha a *1.46 is reduced can be 
exactly as in the laminar case, but the check after the first 
half -step can be eliminated since no stop of the computer oc- 
curred even for Ha2<1.46. 

Reasons for continuation of the calculation may remain 
even after the turbulent separation, but the size of the re- 
maining step must then be readjusted. Often the friction drag, 
after Squire and Young [21], is calculated after a boundary layer 
calculation. This formula is valid only up to the separation 
point, of course. For technical reasons of programming it is 
simpler, however, not to calculate the drag part immediately 
at the separation point but through 

to provide for continuation of the calculation after the sepa- 
ration, in a way that leaves the drag after Squire- Young 
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imcliaiiged. At l^e end ol the entire calculation it can be 
^tained uniformly > 

c« Boundary Layer Shift 

The boundary layer shift is considered in the present 
method simply by using equations (33) and (35), instead of 
i.(16) and (it) for calculations from a certain point on. A 
auddem increase in Hat .at tbe shift point, as was required in 
other methods to soie extent [7], is eliminated. Later on 
the examples will show that in most cases, but not always, a 
strong increase in Has occurs after the shift. It seems sig* 
nificant that this increase is the result of the relations 
used only and must not be introduced artificially. Those cases 
where no rapid increase of Ha 2 occurs appear also of importance 
in practice, which is pointed out during discussion of the ex- 
amples. 

Where the shift takes place is of no special relevance to 
the method discussed. It should be pointed out, however, that 
the most varied shift conditions can be introduced in a simple 
manner. The laminar separation point will be introduced as 
the latest shift point. It has been shown to be a good prac- 
tice to continue the calculation from the laminar separation 
point on principle as for the turbulent case. Additional 
simple shift conditions that are often used are shift when 

U 2 - Ui < 0 (UBl) 


or 


U 2 - Ui - 0 


(UB2) 



That neaas shift at Biniaua pressure. A difference between 
(UBl) and (UB2) exists only when U is constant over a finite 
distance. The transition fron (UBl) to (UB2) pernits calcula- 
tion of the influence, in this case of the advancing nigration 
of the shift point from the end to the beginning of the distance 
with, constant velocity. This influence is ioportant for laninar 
profiles. It detemines the depth of the so-called laainar de- 
pressions in the polar diagram. 

It will be ii^rtant for aany calculations to consider 
the instability of the laainar boundary layer. It depends 
aainly on the fora parameter and on the Re nuaber. The results 
of various stability calculations are plotted in Fig. 8. There /236 
is soae dispersion, depending on the exact shape of the boun- 
dary layer profiles on which they are based. The boundary of 
stability is higher, for instance, according to the mean by 
Ulrich [22], for the Iglisch-profiles of homogeneous bleed-off 
than, according to Pretsch [23], for the Hartree-profiles. A 
certain amount of dispersion is contained in the n thods of 
calculation as well because two different values are given for 
the stability boundary of the so-called Blasius-profile, which 
is mentioned in both categories as a special case. Wieghardt 
[6] carried out a rough interpolation for his calculated values, 
therefore, which has been included in Fig. 8. 

The stability boundary alone does not suffice for predeter- 
mination of the shifting point, however, since a certain dis- 
tance is always required between the onset of instability and 
the shift. The unstable waves must reinforce themselves con- 
siderably at first until their amplitude is sufficient for the 
shift over into the completely irregular turbulent motion. 

During this unstable interval of operation 62 usually increases 
to such an extent that the motion moves, in Fig. 8, from the 
stability boundary away upward, becoming more unstable. For 




Fig. 8 Theoretical Instability 
Boundaries » Experimentally De- 
termined Data for Boundary 
Layer-Shift and Introduced 
Shift Conditions. 


a constant form parameter, 
the Blasius flow for instance, 
only 62 will increase but in 
many cases of pressure in- 
crease a decrease of Hs2 will 
accompany the increase in 62. 
The shift will take place, 
therefore, somewhere above 
the stability boundary in 
Fig. 8 with the degree of 
turbulence of the outside 
flow and the wall roughness 
playing a certain role, of 
course. 

Experiments provide a 
certain indication about the 
practical occurrence of the 
shift for the least distur- 
bance of the flow and the 
wall. Some experimental re- 
sults have been entered in 
Fig. 8 for that reason. 


Schubauer and Skramstad [ 24 ] made their careful experiments 
with the flat plate, i.e., without a change in the form param- 
eter. Granville [ 25 ] investigated the shift for laminar pro- 
files and found a rather clear relation for the difference in 
Re numbers for instability and shift and the mean form param- 
eter during the unstable interval of operation. 




By exchanging the mean form parameter with that at the shift 
point, which is, strictly speaking, only correct for a constant 
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fora parameter, data points from Granville can be plotted in 
Fig, 8 . This was done only for two characteristic values. More 
recent experiments by Raspet [26], for free flight for instance, 
show that the shift occurs even later for smaller Re numbers. 

Two data points from Raspet *s experiment are plotted in Fig. 8 . 
On the other hand, it must be expected that the region of high 
Re numbers and, with it, high form parameter values can only be 
reached through bleed-off. The relation between H $2 and 62 may 
change completely in that case. For example, bleed- off may 
proceed in such a way that the motion takes place always along 
the instability boundary (see examples in section 8 ) or paral- 
lel to it so that there is no move away from the boundary of 
stability during the unstable interval of operation in Fig. 8 . 

It must at least be taken into account that the shift takes 
place close to the boundary there. For that reason the straight 
line plotted in Fig. 8 

In Re 62 > 34.2 H 32 - 46.78 (UBS) 

is recommended as shift condition if one wants to hit the shift 
approximately as it occurs in practice. But if a certain pre- 
vention of the shift is required then the also plotted straight 
line 


In Re 6 a > 34.2 H 32 - 47.81 [U0^ 

is used. If no shift appears in a calculation with this condi- 
tion then one can be reasonably sure that the experiment will 
also remain laminar. 

The choice remains, of course, to introduce different 
shift conditions. Each coidition introduced will normally be 
checked before the calculation of an integration step. An 
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exact check is not necessary for hand calculations until it 
becoaes apparent that the boundary is being approached. 

8 . A Few Bleed-Off Theorems 


The principal boundary layer data can be calculated with 
the methods described, when the normal velocity vo(x) is given 
in addition to the always provided functions Re and U(x). 

Often it is also the other way around, namely vo (x) is to be 
calculated so that the boundary layer will fulfill certain 
conditions} for instance, the prevention of shift of the boun- 
dary layer or of the separation of the turbulent boundary 
layer. 

This problem is also solved easily with the equations de- 
scribed and with the associated methods of calculation and 
programs. Several different aims of bleed-off are jointly at- 
tained when the requirement is to fix the bleed- off velocity 
Vo so that 

H 32 = tK 62 ) * a + b In Refi 2 (51) 

Constant Hs 2 for the laminar and for the turbulent case is 
designated by b=0. But if one chooses a=1.40, b=0.2924, then 
one is always right below the condition for certain shift in 
the laminary case of Fig. 8 and so prevents the boundary 
layer shift. In the turbulent case there are also a few possi- 
bilities with bj^O, which are explained through examples. 

Using ( 11 ) in (51) we get 


''3 = {« I * \ 


(52) 


and through derivation 
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(53) 


g- (. + »(t + 1. «»*)) S+ * vi- f • 

With these two relations the functions 63 > d 6 s/dx and dfia/dx 
can be eliminated from (36) and (37) and one gets 

(4 f (54) 

This makes Vg a function of dU/dx and 62 ; with 6 a one also gets 
H 32 because of (51) and the required boundary layer functions 
can be calculated as before. For continued calculation only 
the change of 62 , given with (36), is required and the differ- 
ential equation (37) is not needed. This exactly corresponds 
to the accepted procedure as given by Wieghardt [ 6 ], for 
instance. 

This procedure does require, however, that the equation 
be already fulfilled at the spot where the calculation with 
the bleed off theorem (51) is to start. But normally* values 
for the pair of functions 62 * H32, are already present from a 
previous calculation at that starting point, which may not be 
compatible with (51) . That difficulty can be eliminated by 
treating (54) and (51) together only as additional equation 
to Vo ( 62 ) and by carrying out the actual boundary layer calcu- 
lation again with the two equations (36) and (37), with no need 
for agreement between the H32 used and 4 >, Because this pro- 
cedure uses the bleed-off theorem only for Vg and not for H32, 
it is questionable whether and how fast H32, which is generated 
in the course of the calculation, adapts itself to the value of 
^ per (51). As the examples will show, that happens very well 
and rather quickly. 

This procedure does require a bit more calculation effort 
than that by Wieghardt, since boundary layer functions for the 
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values 4* aud Has must be calculated. As long as the two values 
are not equal this double calculation cannot be av.. '.ded. When 
calculating by hand equation (37) can be dropped as soon as 
there is close enough proximity between 4* and Hs 2 . For pro- 
gramming this transition is not necessary. This keeps the 
program in its old form, which is beneficial for functional 
clarity. Possible time savings in calculation do not justify 
the effort involved. Besides, the condition vo*0 is possible 
for the recommended method, which is appropriate to many prob- 
lems. According to it bleed -off takes place only when the bleed- 
off theorem really requires it and not when, for instance, 
the boundary layer is brought to the edge of the stability 
boundary artificially, through blow*off,when normally it 
would still be far away perhaps near the stagnation point. 

9. Examples /238 

It should again be pointed out that the unit of length in 
all following examples is that which is also used in the 
Reynolds number Re, according to equation (43). The same 
holds true for the unit of velocity. If one wanted to reintro- 
duce the designations 1 and later for the units, then all 
functions would have to be represented as ratios to these units. 

But that is superfluous. The only thing to be kept in mind is 
that for x»l the unit of length, 1, is meant, from (43), which 
can also be designated as 1; the same is true for all veloci- 
ties . 


a. Parabola- shaped Velocity Distribution for Different 
Step Sice's 

To provide an example for the influence of linear interpo- 
lation on U within one integration step, the curved distribu- 
tion U(x)=x(2-x) for Re*10® with widely varying step sizes was 
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Fig. 9 Influence of Step Size on Pulse Loss Thickness 62 

calculated. Consecutively, Ax»l, 0.5, 0.2,0.1,0.02 and 0.01 was 
chosen. Results for the four biggest step sizes are shown in 
Fig. 9. The comparison was demonstrated only on 62 . The errors 
in H 32 are smaller. The data points for U used and those calcu- 
lated for 62 were always joined by straight lines to character- /239 
ize the step size. The effect of decreasing Ax to below the 
value of 0.1 can no longer be shown graphically. For that 
reason the results were put together in Table II for the two 
points x“l and x* 1.6 with their errors. 
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The step size 0.1 elreedy sho«s food results when a total 
of only sixteen steps is counted in the interval. But gross 
errors mist be expected when the velocity increase is linearized 
too roughly and the stagnation point solution too far extended, 
which often happens. 


Table II 


Pulse Loss Thickness for lf*x(2'X), Re«10* for x«l and 
x«1.6 as Functions of Step Size 
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b. Howarth-Flow with Turbulent Continuation 


The velocity distribution U(x)*l was treated rigorously by 
Howarth [27] and recalculated accurately by Leigh [28]. This 
calculation often served for conparisons with the results fron 
approxination nethods (halt [5], for instance). In particular 
the point of laninar separation is very critical and could 
heretofore not be calculated with satisfactory accuracy by ap- 
proxiaatinn nethods. The equations and functions now available 
produce an error that is smaller by an order of nagnitude. A 
graphic conparison with the exact solution is no longer pos- 
sible because the lines do not provide sufficient difference 
for it. Nunerical values for the separation point are shown 
for various nethods in Table III. It cones as surprise that 
the present procedure offers nuch better results than that by 
Walz even though the basic equations and the boundary line 
profiles on which functions H 32 > e* and 0* are based are the 
sane ones for the laninar case. As mentioned before, it was 
found that the functions used by Walz had been incorrectly 
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Table HI 


Separation Point for Howarth-Flow in Various Directions 


Procedure 

X 

(Separation) 

Error 

% 

Howarth* Leigh 

0.1198 

0.0 

Pohlhausen P4 

0.16 

33.0 

Pohlhausen with Hartree profiles 

0.103 

14.2 

Walz with Pohlhausen P4 

0.125 

4.2 

Walz with Hartree profiles 

0.114 

5.0 

Wieghardt 2 parameter method 

0.116 

3.3 

Gortler (difference method) 

0.1183 

1.2 

Prel. method with Walz-Hartree function 
Prel. method with new Hartree functions 

0.1137 

5.0 

AX«0.01 

0.1202 

0.4 

Ax»0.002 

0.1199 

0.1 


calculated. With the interrelations for H 32 , e* and D* used 
by Wal '7 the sane results are achieved here as there, with the 
present fomulas. Further inprovenent can only be attributed 
to accurate evaluation of the Hartree profiles. The actual 
stateneat by Walz is significantly better, therefore, as can 
be concluded fron the results up to now. 

Fig. 10 also gives the curve of the fom parameter H 32 for 
various Re numbers. The laminar separation point was here 
assumed as shift point. In the laminar part Hsj is independent 
of Re, which can already be deduced from the basic equations, 
but that is not true in the turbulent part. In Fig. 10 it can 
be clearly recognized how H 32 increases rapidly with the higher 
Re numbers after the shift, at approximately the rate at which 
Truckenbrodt [7] introduced it artifically. But for smaller 
Re numbers the jump does not occur. Here Ha 2 still moves for 
awhile after the shift in the region H32<1.58, where separation 
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is already possible. This result agrees with the observation 
that local separation is often found (29) before the boundary 
shift, for smaller Re numbers. According to Fig. 10, the 
reason for this separation bubble must not only be sought in 
the insufficiency of excitement of the unstable disturbances, 
but also in the too slow formation of the turbulent boundary 
layer and with it the removal of flow reversal danger at the 
wall for this range of Re number and form parameter. 

Additional calculations with uninfluenced turbulent boun- 
dary layers have already been carried out in large numbers, 
also comparisons with experiments. Two of these comparisons 
have been published previously (14), There it is shown, in 
particular that the increase of <52, according to the much 
applied method by Truckenbrodt for the area of high pressure 
increase, seems too small while better values are achieved 
with the equations presented here. One considerable improve- 
ment of Truckenbrodt ' s method ,mer t ioned by Schol: (10), does 
correct the shortcoming for the pressure increase in the turbu- 
lent region but generalization of this improvement for the 
laminar case and for bleed -nff has not been tried so far. 

c. Laminar Plate Flow with Homogeneous Bleed-Off and 

A further classical example, for which an exact solution 
is available, is the flat plate 0=1, with constant bleed-off 
Vo. The recalculation of this case, which is shown in Fig. 

11, also provides very accurate results. Iglisch's ( .'0 1 solu- 
tion agrees in ^2 almost completely with the present one, 
there being only insignificant differences for Il32- To achieve 
such good agreement the step size of the first step Ax, which 
is calculated without bleed -off here, must be kept very small 
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tter Hs 2 for Howard 
‘ent Re numbers. 






since the influence of bleed-off is particularly significant 
there. 


For conparison, the lengths that were nade diaensionless 
by Iglisch are recalculated. In the function used there 




' I 


only the Re nuaber need to be introduced «per (43), along with 
the units 1 and U^; the interrelation between the values chosen 
here for v« and Re, and C and x, can then easily be established. 
The sane holds for function 

f I/. I r • 

The sane case of plate flow but with positive v« , i.e., 
blow off, is shown also in Fig. 11. As in all approximation 
methods, laminar separation occurs for x»0.4, which cannot 
occur for rigorous treatment. That is the case mentioned in 
section 5(b) in which the present method does not work too well. 
In practice it is unimportant since the blow-off profiles are 
anyway very unstable and boundary layer shift cannot be avoided 
here. Still, it may serve as an example of how little signifi- 
cance errors found in special cases have for the entire method. 
As far as the laminar part is concerned in the present method, 
it must be concluded from results up til now that it offers 
very good results when pressure increase without bleed-off or 
bleed-off to high values of H $2 is involved. This exactly 
corresponds to the choice of functions Hia, e* and D*, made 
in section 3. For more detailed investigations of errors we 
refer to Nickel (31). 
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d. Plate Bleed-Off with Constant Form Parameter 


In Fig. 12 the distributions 62 (x) and Vo (x) for the flat 
plate with constant form parameter, which have been obtained 
with the present method, are compared with the exact solutions. 

Again very good results are generated for H32>1. 57258. That is 
only natural because the form parameter interrelation was taken /242 
from the corresponding exact solutions. On the other hand for 
blow-off, i.e.j for H32 values below 1.57258, results would be 
somewhat worse because for that region other in*‘errelations for 
form parameters were chosen. Since blow-off is not allowed in 
the programs used here, for reasons already mentioned in sec- 
tion 8 , blow-off cases of the "similar" plate flow (anyway 
again largely unstable) were not calculated. 



Fig. 12 Recalculation of Bleed Off Velocity vo for Similar 
Bleed-Off Boundary Layers at U»const., Re*10®. 
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e. Bleed-Off of Laminar Boundary Layer for Strong Pres- 
sure~Tncrease 


An important aiea of application for bleed -off is the 
prevention of boundary separation, for instance on air foils 
with high lift coefficients. Large negative pressures and, 
subsequently, long and strong pressure increases occur on the 
topside. To get an overall picture for the bleed-off required 
in such cases a large number of calculations were carried out 
for a velocity distribution composed of linear sections, shown 
in Fig. 13. The pressure increase in the distribution on which 
it is based goes from U=4 to U=0.4, thus attaining a value that 
could never be overcome without bleed-off. This order of 
magnitude approximately corresponds to what occurs on air foils 
with very high lift coefficients (between 3 and 5). 


ie% 



Fig. 13 Quantity Coefficient Co and Plus Loss Thickness 62 
for Strong Pressure Increase and Form Parameter of the Laminar 
Boundary Layer, H 32 , Being Constant; R = 2.10®. 


In Fig. 13, as well as in the following Figs. 14-16, the 
function 




X 

» 

Vodx 


(55) 
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is used throughout in place of vo, in addition to and 
This function offers a better overall picture than vo itself 
since the total amount of bleed-off is more important. At the 
trailing edge Cq agrees with the often used quantity coeffi- 
cient Cq, before that c^q provides the amount of bleed -off 
needed up to a point x while -vo is proportional to the slope 
of curve Cq. 

A simple method for preventing separation is holding the 
form parameter constant through appropriate bleed-off. This 
bleed -off calculation can be carried out for laminar or turbu- 
lent cases by means of ( 54 ), setting b* 0 . At first the laminar 
boundary layer was calculated for various form factors, without 
consideration of instability and shift. The resulting values 
for 62 and Cq are shown in Fig. 13 . At first glance the re- 
sult is surprising. One suspects at first that bleed -off to 
high H32 values should be much stronger and give correspond- 
ingly smaller values for 62. The calculation offers for x»l 
values of Cq that differ very little, while 62 is high for 
high H32 values and low for low ones, the exact opposite of 
what is expected offhand. Only at the start of the calculation 
does the picture correspond to expectation. As soon as this 
condition has become pronounced a smaller 62 occurs because of 
( 54 ), in consequence (and because of the larger e*) a larger 
dfi2/dx with the increase in 62 predominating, so that the pic- 
ture is quite reversed in the end. The reversal is easily 
recognized when ( 54 ) is entered into ( 36 ), b »0 and H3 2*<t». 

This results in 

2/>»— f* /•|_Wi,— 1\. 1 dl 

dx n„~'ipv dx’ ( 56 ) 

The first term, which is not dependent on 62dU/dx, actually 
gives a smaller dfi2/dx with increasing H32 since the numerator 


/243 


47 



decreases and the denominator increases, according to Fig. 2. 

The second term shows the opposite tendency. Since 62 is in 
the denominator of the first term, but in the numerator of the 
second one, the character of Fig. 13 is easily understood. 

The second term predominates as soon as 62 has become larger. 
Since 62 also goes into Vg the lines for Cq diverge later, but 
much less so than those for 52 - 

During prevention of the separation the required amount 
of bleed -off is predominant while 62 only goes into the resis- 
tance ianyvay quite high for high life coefficients). Ac- 
cording to Fig. 13, the value of the form parameter, up to which 
bleed -off takes place, does not have any big influence on the 
amount of bleed -off. 

That result is important because it shows that stability 
of the laminar boundary layer, which increases strongly with 
Ha 2 , can be obtained even for strong increases of pressure 
without a very large amount of bleed -off. This also confirms 
another calculation where two pairs a,b, are entered into (54). 
For case (A) a*1.40, b*0. 02924, was chosen so that the region 
below the strict shift condition in Fig. 8 not be exceeded, 
while case (B):a=1.37, b=0. 02924, permits the region to 
the upper straight line in Fig. 8 , along which shift may al- 
ready be expected. The results shown in Fig. 14 again demon- 
strate very small differences for Cq, but bigger ones for H 32 
and 6 2 . 

Particular attention is called to H 32 having a maximum. 

That points to a maximum of Refi 2 , which is a function of the 
strong increase of U. In case (B) even 62 has a maximum. 

Since -Vo also becomes very big, along with function (l/U)dU/dx 
at around x=l, 62 does not continue to increase here any longer. 
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Fig. 14 Quantity Coefficient cq, Form Parameter H 32 and 
Pulse Loss Thickness 62 , for Strong Pressure Increase, Re = 
2.10®- -(A) Bleedoff for Certain Prevention of Shift, (B) Pre- 
vention of Shift no Longer Certain. 

By putting (54) into (36) a closer picture can be obtained 
here too. 
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f. Bleed-Off of the Turbulent Boundary Layer at Strong 
Pressure Increase 


Since it is often not easy to maintain the boundary 
layer laminar with certainty, various attempts have been made 
recently to bleed-off the boundary layer so that no separation 
could occur [ 2 ]. To get some indication about how correct 
and rational such bleed-off is, some examples of turbulent 
boundary layers were calculated also for the velocity dis- 
tribution of the previous examples. 

Here, too, it is again easiest to start calculating with 
constant form parameter. The resulting values for Cq are 
shown in Fig. IS. Bleed-off starts a little later than for 
the laminar cases because there is a wait for the shift. 
Compared to the laminar cases the picture shows a complete 
shift. The Cq lines show surprising overlapping. For smaller 
H32 values they seem to be more favorable. Later on the per- 
mitted increase of 62 during the initially small Vq takes its 
toll with bleed-off amounts that can hardly be realized. 

The high values for H32 on the other hand are very high at 
the start with the required bleed-off. Even though bleed-off 
with a value of H 32 * 1.95 is most favorable in the end, the 
case appears unfavorable when compared to laminarity. The 
impression is gained therefore, at first, as if turbulent 
bleed-off is much less favorable. A series of additional 
calculation of samples shows, however, that the difference 
need not be so great. Apparently, bleed-off must be carried 
out with variable form paramter so that at first the values 
for H32 are small, increasing later on. For that reason ( 54 ) 
was used for the calculation and constants a,b set so that 
for, 
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In » 10(Re^^*2000) Haa • 1.95, 

and for 

In RCg^ » 6(Reg^=400) H 32 - 1.95 + k. 

That results in 

a - 1.95 + 2.5 k b » -0.25 k. ( 57 ) 



Fig. 15 Quantity Coefficient cq and Pulse Loss Thickness 62 
for Strong Pressure Increase ana Constant Form Parameter H 32 
of the Turbulent Boundary Layer, Re » 2.10®. 

Since for higher values of k the form parameter attains 
high values close to 2, the calculation was carried out with 
a,b, calculated from (57) but only to H32*1.99; above that 
a»1.99, b»0,was used. 

The results are shown in Fig. 16. It appears that values 
for c^Q have decreased significantly. It should be remembered 
that these calculations stop at values for H 32 that are far 
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away from separation. Much more bleed-off than was necessary 
for prevention of separation took place then. By takinj^ that 
into consideration and bleeding less off the trailing edge Cq 
will perhaps not decre..se exactly down to laminar v'alues but 
it will be of the same order of magnitude. 

The most favorable combination in Fig. 16, a=1.35, b=C.06, 
also turned out to be favorable for other cases where velocity 
is not linear. For the choice of a,b, in the turbulent case 
the total pressure increase is almost the only decisive indi- 
cator. 

Values of Cq obtained so far are of the same magnitude as 
those obtained experimentally by Raspet [2]. Unfortunately, 
no straightforward comparisons are possible since- no experi- 
ments have been published so far for which distributions of 
pressure and bleed-off have beer supplied. Another difficulty 
should be pointed out The turbulent bleed-off calculations 
may in part lead to combinations of H 32 and 62 where, according 
to estimates from section 4, the boundary layer thickness ex- 
ceeds by far anything that is expected from the derivation of 
the boundary layer equations. Raspet and Cornish have pointed 
(unpubl _.hed) to experimental results that are not in agreement 
with the boundary layer equations. Understandably ,62 increases 
less in the experiment than in theory. On the one hand the 
theoretically available contribution F to 62 , according to (30)i 
cannot be determined experimentally from regions where the dif- 
ference between u and U is very small, since potential flow can 
already be detected in such regions. 

On the other hand the contributions of 62 that exceed a 
permissible boundary layer thickness reach into regions of 
potential flow where dU/dx is smaller than at the wall. It is 
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to be expected that the i»'crease of 62 , according to the pulse 
theorem, will not be reached, i.e., that the theoretical values 
for 62 and with it for Vo, are too large. That wjII keep me 
on the safe side for most problem situations. 



Fig. 16 Quantity Coefficient cq. Pulse Loss Thickness 62 , and 
Form Parameter H 32 for Strong Pressure Increase and Bleed-Off 
of the Boundary Layer per (5^) and (56), Re=2.10®. Dashed line: 
Cq of case (A), Fig. 14. 

Turbulent bleed-off is also to be considered as definitely 
positive for the high-lift case in general. The q-j ,tion which 
of the two types of bleed-off is preferable must be decided by 
practical considerations. Such a significant practical consid- 
eration would be that of safety. 

Once a separation has taken place it can only be removed 
by a basic change in the pressure distribution, in all cases 
of bled off boundary layers. In the polar diagrams of air 
foils that produce a marked hysteresis both in the laminar and 
in the turbulent case. But judging from the line for the 
laminar Cq in Fig. 16 tnat requiies a little less bleed-off. 

But boundary layer shift alone can also lead to separation. 
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Since the high-lift case occurs priaarily during takeoff and 
landing, where gnats, bugs and birds can force a boundary layer 
shift, it seens nor^^ppropriate here to apply turbulent bleed - 
ofC ind even te aoniaor constantly whether laninar flow was 
not generate tbeough oversight. In theory that can hardly 
happen siac« bleed -off cones too late for it. In practice 
bleed off i;> .ore or less discontinuous so that the slight 
differences bethreen laainar and turbulent siphoning are hard 
to verify at the start. It will probably be required to take 
precautionary neasures not to work in a region where pure 
boundary layer shift alone can already cause separation and 
cr»sh. It is possible that this nay be the reason for two 
serious crashes that are so far unexplained. 
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